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We investigate the structural conditions for metamagnetism in MnP and related materials using 
Density Functional Theory. A magnetic stability plot is constructed taking into account the two 
shortest Mn-Mn distances. We find that a particular Mn-Mn separation plays the dominant role in 
determining the change from antiferromagnetic to ferromagnetic order in such systems. We establish 
a good correlation between our calculations and structural and magnetic data from the literature. 
Based on our approach it should be possible to find new Mn-containing alloys that possess field- 
induced metamagnetism and associated magneto caloric effects. 

PACS numbers: 75.30. Sg,75.30.Kz, 75.80.+q, 75.30.Et 



INTRODUCTION 



The magnetocaloric effect (MCE) is the change of tem- 
perature of a material in a changing magnetic field, and 
was first observed in iron by Emil Warburg. It has until 
now been exploited in paramagnetic salts as a means of 
cooling below the temperature of liquid Helium. How- 
ever, the high global warming potential of conventional 
HFC refrigerants^ and the large fraction of energy used 
for cooling are concerns that currently fuel interest in 
such "magnetic cooling" as a high-efficiency solid-state 
method of refrigeration close to room temperature. In 
recent years much of the search for viable room tem- 
perature magnetic refrigerants has focused on low cost 
transition metal-rich compounds that possess first or- 
der magnetoelastic or magnetostructural phase transi- 
tions that can be driven by applied fields of the order of 
1 Tesla. Most of these materials are Mn-containing fer- 
romagnets such as Ni2MnG8P, Ri_xMnxMnO#l, MuAP 
and MnFe(As,P)^. The observation of MCEs during 
order-order transitions (in metamagnets) is more rare 
and the main examples of l^ir^ effects near room tem- 
perature are in FeRh^^, CoMnSi^and Ni2MnSrP. In some 
ways metamagnets offer a natural route to the desirably 
elevated MCEs associated with a first order, rather than 
continuous magnetic transition. However, there are few 
metamagnets examined for their MCE that have not al- 
ready been synthesised for other reasons. In this pa- 
per, we focus on MnP-type {Pnma space group) binary 
and ternary orthorhombic materials and in particular on 
what structural factors support antiferromagnetic order 
and field-induced metamagnetism. Our motivitation is 
to allow prediction of new materials that may have large 
MCEs. 

The appearance of rich magnetic features with pecu- 
liar non-collinear (NC) magnetism such as fan, helical or 
cycloidal spin structures observed in MnP-type binary 
and the related XMnZ (TiNiSi-type) ternary alloys has 
engaged the interest of both theoreticians and experi- 
mentalists for several decades^^"^^". In what we term the 
"prototype" MnP binary alloy, the appearance of the so- 
called "double-helix" spiral state has been the subject of 
several studies. This is a state in which two pairs of 



Mn atoms in each cell take on the same helical wavevec- 
tor, but with a fixed phase difference between the pairs. 
Gribanov and Zavadskii^^ ascribed the occurence of the 
magnetic double helix in MnP-type phases to the exis- 
tence of indirect double exchange coupling. Based on an 
analysis of the energy stability of helimagnetic spin ar- 
rangements, Fjellvag et al. discussed the possibility of 
a metamagnetic transition by purely carrier localization 
with no changes in structural parameters. The compli- 
cated crystal structure of MnP is, however, a serious ob- 
stacle to a complete analysis of magnetic exchange inte- 
grals based on Goodenough's double exchang^^. Kallel 
et al.^^ concluded that the presence of high ratios of an- 
tisymmetric to symmetric exchange is necessary to ob- 
serve helical structures. They constructed a model with 
exchange integrals from only the nearest and second- 
nearest Mn-separations taken into account. Dobrzynski 
et. al.t^ argued that it is necessary to consider the rela- 
tion of exchange interactions with up to the 7th nearest 
neighbours in order to account for the magnetic structure 
of MnP. In TiNiSi-type ternaries, a similar complexity 
emerges. Studies of Coa^Nii.xMnGe have similarly fore- 
cast the need to consider up to 8 exchange interactions, 
especially when a non-zero moment is present on the Co 
site^^ 

Despite extensive studies of the formation of non- 
collinear magnetic phases in MnP-based materials, no 
simple, guide-like explanations have been deduced that 
may allow the design of new materials with a field- 
induced metamagnetic transition near to room temper- 
ature. In this article, we compare our theoretical cal- 
culations, which use density functional theory (DFT), 
with experimental findings to depict the most important 
structural parameters for practical material design. Our 
observations of high MCE and giant magneto-elastic cou- 
pling in CoMnSi^ and the rich variety of materials in 
the literature that take on the MnP (or TiNiSi) struc- 
ture^^ motivate our study of metmagnetism in the MnP 
class. 

The paper is organised as follows: in sec. [Hj the un- 
derlying Pnma crystal structure is introduced together 
with the details of our computational method. Sec. [HI 
describes the results of ab-initio electronic calculations 
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Figure 1: (X)MnP-type {Pnma, 62) orthorhombic structure 
(the diagram was visualized using VESTA^. The domi- 
nant di and d2 Mn-interatomic distances are constructed with 
dashed and sohd hues, respectively. 



based on DFT theory. We will show that the ferromag- 
netic structure of the prototype MnP binary alloy would 
undergo a change in its magnetic state as a result of isodi- 
rectional lattice expansion. This change, however can 
also be realised by chemical pressure through alloying as 
is demonstrated in Sec. |IV[ Examples from the literature 
for both 3d transition metal addition as well as p-block 
partial replacement in MnP will be discussed. In Sec. [V| 
we extend this analysis to other Mn-based ternary (TiN- 
iSi structure) compositions that crystallize in the Prima 
space group. Finally, we summarise the practical factors 
that can scale the crystal structure into the regime where 
a change of magnetic state can be expected. 



II. CRYSTAL STRUCTURE AND 
COMPUTATIONAL DETAILS 

A. Crystal Structure 

The orthorhombic crystal structure of MnP (Pnma, 
62) is shown in Figjl] with unit cell dimensions 
a=5.268 A, 6=3.172 A and c=5.918 A at room tem- 
perature. Both the Mn and the P element occupy 
the 4c (x, |,z) crystallographic positions with xmu = 
0.0049(2)^Mn = 0.1965(2) and xp = 0.1878(5), zp = 
0.5686(5p^. The orthorhombic structure can be regarded 
as a distortion from the higher symmetry hexagonal NiAs 
(P63/mmc, 194) type structure with the atomic posi- 
tions of Ni (4c) shifted from xmu = 0, z^i = ^ to the 
ZMn positions and the As (4c) atoms moved from 
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Table I: Possible collinear magnetic configurations of 
Mni_4in a single unit cell of MnP and their relation to di^d2 
inter-Mn distances. 

xni = \, ZNi = j2 ^P' Furthermore the two 

structures can be related as follows: , bortho = cihex and 

^ortho 

= ^/3 X ahex- One direct extension of the MnP- 
type binary structure is the TiNiSi-type ternary, with the 
occupation of an additional atomic species at a general 
4c crystallographic position as shown in gray in Fig. [T] 
Therefore in the case of a general XMnP alloy, the num- 
ber of minimal basis atoms is 4 for each species in the unit 
cell and assigned as Mni_4 and P1-4 for a binary with 
an additional Xi_4 element in case of a ternary alloy in 
Fig. [1] In addition, the two smallest Mn-Mn distances, 
the focus of our attention in later sections, are denoted 
as di and ^2- 



B. Computational Details 

The electronic structure calculations were carried out 
using the VASP code, based on density-functional theory 
(DFT)^^. Site-based magnetic moments were calculated 
using the Voskown analysis within the general gradient 
approximation (GGA) scheme. We used the minimal, 8 
atom basis supercell for the binary MnP alloy with 4 Mn 
atoms and 4 P atoms on the 4c crystallographic posi- 
tions, respectively to calculate the total energies of the 
ferro-, antiferro-, and non-magnetic states. The accu- 
rate DoS with high density of k-grid of 21, 19, 21 points 
along the a, b and c-axes were used, respectively. The 
spin-orbit interaction was turned off during the calcu- 
lations and only collinear magnetic configurations were 
considered. The energy convergence criterion was set to 
10~^ eV during the energy minimization. Where applied, 
the expansion and compression of lattice parameters were 
realized isotropically. 

III. RESULTS: MNP 

As outlined in Section [T[ MnP exhibits complex non- 
collinear ferromagnetic features including low temper- 
ature screw and conical magnetic structures, depend- 
ing upon the axis of magnetisation. These transform 
into collinear ferromagnetism at higher applied fields and 
temper at ure4^. The collinear ferromagnetic state is sta- 
ble above T > 47 K and it has a Curie temperature 
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Figure 2: (Color online) Total energy of non-magnetic (NM), 
ferromagnetic (FM) and antiferromagnetic (AFM) solutions 
as a function of isotropic lattice expansion and compression 
in MnP. 



of Tc = 290 kP. Such peculiar non-collinear mag- 
netism was explained using spin wave theor}^^'^ and was 
more recently ascribed to the influence of magnetocrys- 
talline anisotrop}^^^. There have been several theoreti- 
cal studies of both the magnetism and crystal symmetry 
of MnP. The first theoretical work based on augmented 
plane wave (APW) theory reported ferromagnetic ground 
state with moments of 1.2/iB/f.u.'^, in good agreement 
with the experimental value of 1.29/i^ . Focussing on 
the possibility of NiAs-MnP structural distortion, Tremel 
et. al. showed the importance of direct and indirect 
metal-metal bonding due to a second-order Jahn- Teller 
type distortion within a tight-binding (TB) approach^^. 
They found that the stabilization of the lower symmetry 
MnP phase over the hexagonal one is due to an increased 
metal-metal interaction, which results in a lower density 
of states at the Fermi level. Spin polarized calculations 
on the stability of the orthorhombic phase over the hy- 
pothetical NiAs and zincblende-type ones in MnP was 
recently also investigatecF^^ and found good agreement 
with experimental lattice values. 

In this article, we fix the lattice type as Pnma as 
found experimentally and concentrate on the electronic 
and magnetic properties of the prototype MnP binary al- 
loy. We calculate the effect of isotropic lattice expansion 
and compression on non-magnetic (NM), ferromagnetic 
(FM) and antiferromagnetic (AFM) solutions in order to 
find the critical lattice parameters where the crossover 
from one magnetic state into another can occur. For this 
study a single unit cell consisting of 8 atoms (4 Mn and 4 
P) is used, which allows three different collinear antiferro- 
magnetic configurations (AFMl, AFM2 and AFM3) and 
a collinear ferromagnetic (FM) one to be constructed, as 
show in Table m 

Fig. [2] shows the calculated total energies of the con- 
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Figure 3: (Color online) Difference between the total en- 
ergies of antiferromagnetic (AFM) and ferromagnetic (FM) 
solutions in MnP as a function of di and (i2, the two short- 
est Mn-Mn separations. AFM configurations become stable 
where AEafm-fm < 0. The vertical dashed line represents 
the experimental (strain-free, e — 0) Mn-separation values. 



verged solutions as a function of isotropic lattice expan- 
sion and shrinkage. As is reflected from the figure, the 
collinear ferromagnetic state is predicted to be energeti- 
cally most favorable in accordance with experimental ob- 
servations (above 47 K) in the range —0.06 ^ s ^ 0.05. 
Under severe compression {e < —0.06), where the clos- 
est Mn-Mn atomic distance is d2 ^ 2.55 A, no spon- 
taneous magnetisation is stable in the broad overlapping 
3d bands and the calculations find convergence instead in 
non-magnetic solutions. On the other hand, with lattice 
expansion, antiferromagnetic groundstates can become 
more stable than the ferromagnetic one, as indicated by 
the crossover of the total energy line of AFMl around 
e ^ 0.05 in Fig. [2j Further lattice expansion leads to the 
AFM3 configuration becoming the most stable one above 
e>0.n{d2< 3.00 A). Notably, AFM2 is never favored. 

The relative stability of the different magnetic ground 
states becomes more evident if the differences between 
the energies, Etoi of AFM and FM states are compared 
{AEtoi = Eafm — Efm)' By this comparison, a non- 
FM state becomes most favourable when it has the most 
negative value of AEtoi- In Fig-[3) the regions of stability 
of the different magnetic states are shown as a function 
of Mn-Mn distance. The FM state is most stable for in- 
termediate deformations whilst the non-magnetic mode 
and AFM modes become favored at high compression and 
at high strain respectively. Considerable magnetic mo- 
ment was found only on the Mn site (Mmu) as a result 
of exchange splitting of the d-orbitals with relatively in- 
significant magnetic moments, Mp < O.l/i^, on the sites 
of P atoms. Figure [4] shows the magnetic moment of Mn 
atoms in the different magnetic states as a function of 
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Figure 4: (Color online) Mn magnetic moment vs. strain for 
different magnetic states in MnP (see Table 



lattice size. In the compressed lattice region (left part of 
Figs. [3] and [4| tfie close proximity of Mn atoms causes 
a strong overlap of their (i-orbitals, resulting in broad 
d-d hybrid bands that cannot support spontaneous mag- 
netisation. The magnetisation that grows with lattice 
expansion can be interpreted as a result of smaller d-d 
overlapping that makes the d-orbitals of Mn atoms more 
localized. This G?-orbital localization then enhances the 
exchange splitting of the Mn d states reflected in the 
increased magnetic moments. 



In a compressed unit cell the p-d hybridization of P 
and Mn atoms can also lower the d-d exchange interac- 
tion, reducing the magnetic moments. The enhanced ex- 
change splitting with increased cell size is demonstrated 
in Fig. ^jp (bottom), which shows the total density of 
states (DOS) of the ferromagnetic MnP. In the equilib- 
rium volume {s = 0), the main peaks of Mn c^-band 
in the majority and minority states are separated by 
0.9 eV, which increases to 1.5 eV for the lattice with 
e = 0.055 lattice expansion. Simultaneously, there is 
an increase of the total density of states at the Fermi 
level, NrotiEw) = N^{Ew) + N^{Ey) in the FM state 
from NToti^w) = 5.3 to 5.9 states/eV/f.u. However, a 
lowering of the DOS of the AFMl configuration from 
NTotiEw) = 5.1 to 1.8 states/eV/f.u. can be seen, which 
suggests that the AFM phase is energetically more favor- 
able Fig. |5|l (top). 

In the next sections, we will show that our stability 
plot (Fig. [3|, calculated for the binary MnP alloy can 
be used to predict the stability of AFM or FM states in 
other (X)MnZ type binary and ternary alloys with the 
same (Pnma, 62) structure. 



> 




CO 8 

I ' 

w 4 
o 2 

-8 
-10 



~\ AFMl ' 
^' i I.'"', 

A p I > ; 


1 if 

\ il 






k 


b) 'Z:i^^ MM, 


1 


V wy(\ 




M II \ 



-^E-Erev| 



Figure 5: (Color online) Density of states for binary MnP al- 
loy. The total density of states at the Fermi level (NrotiEw)) 
decreases for antiferromagnetic (AFMl) configurations (a) 
with lattice expansion from £ = (solid line) to £ = 0.055 
(dashed line), while an opposite tendency is observed in (b) 
for FM alignment . 



IV. IMPLICATIONS FOR XMnP AND Mn(P, Z) 

A. (Fe, Co)MnP 

A direct extension of the binary MnP alloy is the 
ternary FeMnP. The Fe atoms can occupy another 4c site 
without changing the orthorhombic crystal symmetry. 
Interestingly, the addition of iron atoms to the ferromag- 
netic MnP results in a commensurate non-collinear AFM 
structure with Tn = 340 K. The magnetic unit cell is 
twice as large as the structural one, being doubled along 
the c-axis with magnetic moments of Mmu = and 
= 0.5/iB lying in the ab-plane^^. The antiferromag- 
netism of the system can be explained by means of the 
stability plot drawn for MnP in Fig. [3| The added Fe 
atoms cause an expansion in lattice dimensions, leading 
to an increase in both di and 6^2 inter-Mn separations into 
the regime where AFM coupling is stable. The resulting 
di = 3.05 A distance is in the proximity of the region 
where both AFMl and AFM3 collinear configurations 
are close to the lowest point in total energy, which indi- 
cates the trend to a complex non-collinear AFM ground 
state. 

The dominant effect of the nearest Mn-Mn distances in 
determining the magnetic nature of these alloys becomes 
more pronounced with the investigation of the pseudo- 
ternary (Co,Fe)MnP composition. (Fei_xCoaj)MnP alloys 
exhibit complex magnetic properties ranging from non- 
collinear AFM for X = through temperature-induced 
metamagnetism for 0.5 < x < 0.8 to collinear FM struc- 
ture above x > 0.^^. Extensive theoretical work based 

I I 

on KKR-CPA calculations by Zach et. al.^ tentatively 
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explained the NTotiE^) to be partially accountable for 
the metamagnetic transition. Their investigation on var- 
ious collinear AFM configurations found a formation of 
a 'pseudogap' between the conduction and valence elec- 
tron bands with Co addition, which results in a lowering 
of NToti^F) similar to that in the strained MnP alloy of 
the previous Section (Fig. [5|. 

The same authors have also tried to associate the 
observed magnetic states to changes in crystal struc- 
ture. They found the measured Mn-(Fe,Co) and (Fe,Co)- 
(Fe,Co) inter-atomic distances to be fairly independent 
of temperature. The metamagnetic transition is ob- 
served in the region where 2.90 A < di < 2.95 A and 
3.15 A < d2 ^ 3.16 A with a change much more visi- 
ble in the di distance as Fe is partially replaced by Co 
when compared to the change in ^2 • In Fig. [6j the inter- 
atomic Mn-Mn distances {di^d2) in the pseudo-ternary 
(Fei_xCoa;)MnP alloy are plotted as a function of effec- 
tive atomic radius. Here the effective atomic radius is 
defined as Reff = x x Rco + (1 — x) x Rfc^ where x 
is the concentration, Rco is the atomic radius of the 
Co atom (152 pm) and Rpe is the atomic radius of Fe 
atom (156 pm), respectively. The metamagnetic transi- 
tion temperature (T^) shows a strong dependence on Co 
concentration. It has a value of Tt ~ 300 K at x ~ 0.5 and 
is completely suppressed at x ~ 0.&^. Comparison of the 
interatomic distances with those seen in the MnP stabil- 
ity plot of Fig. |3] shows a remarkably good agreement 
between the theoretically expected and the experimen- 
tally found regime of AFM-FM metamagnetism. As the 
di interatomic Mn distance is increased above '^2. 95 A 
the overall structure is shifted into the region where an- 
tiferromagnetism is the dominant state. This change of 
magnetism can also be achieved by partial replacement of 
the p-block element as demonstrated in the next section. 



B. MnAsi-xPx 

The ternary MnAsi_xPa; alloy is another example of 
metamagnetism and has been extensively studied by Fjel- 
Ivag et. al. in the < x < 0.20 composition range^^. 
Its rich magnetic and structural phase diagram includes 
the occurrence of a magneto-structural transition be- 
tween hexagonal FM and paramagnetic (PM) orthorhom- 
bic phases around room temperature for x < 0.025. That 
transition, in the related material in Mui.xFe^^As, re- 
cently attracted much attention due to the claims, now 
disputed, of an anomalously high MCE However, the 
interest in the origin of the crystal symmetry change is 
long-standing. Initially, Goodenough and co-workers as- 
sociated the transition of electron states of Mn atoms 
from high-spin to low-spin to the change of lattice struc- 
ture from hexagonal to orthorhombic in Mn(P,Asj'^. Re- 
cently, Rungger and Sanvito using detailed ab-initio elec- 
tronic structure analysis found that if antiferromagnetic 
alignment in the basal plane of the hexagonal structure 
is imposed, the orthorhombic structure becomes more 
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Figure 6: (Color online) Experimentally observed Mn-Mn 
distances {d\,d2) in (Fei-xCocc)MnP pseudo-ternary alloy at 
100 K (aftei*^*^). Antiferromagnetism becomes stable as d\ 
increases above '-^2.95A, in line with our calculations on MnP 
in Fig. [3] Temperature-induced metamagnetism is seen for 
0.5 < X < 0.8. 
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distances (di, ^2) in ternary MnAso.88Po.12 pseudo-ternary al- 
loy as found experimentally*^. In accordance with Fig.jsjthe 
AFM state is stabilised when d\ is greater than ^^2. 95-2. 98 A. 



stabldSH. 

As with MnP case, we focus on magnetism within 
the orthogonal, Pnma structure. Compositions around 
X ^ 0.10 show complex magnetic behaviour with domi- 
nant Ha-typ^^ helicoidal AFM spin arrangement at low 
temperature that reappears at high temperatures with a 
transient collinear FM configuration in between. Fjellvag 
and co-workers recognised the major importance of the 
shortest Mn-Mn separations and used neutron diffrac- 
tion to map them as a function of temperature. Fig. [7| 
reproduces the temperature dependence of di and d2 
that they found in MnAso.88Po.i2- Interestingly, while 
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Figure 8: (Color online) Mn-Mn separation (di, ^2) observed 
in XMnSi ternary alloys, as a function of effective X atomic 
radius (low temperature data)^^ 36j|i39j j^Q^gj^g^j lines are fitted 
linearly to data and are only a guide to the eye. 



the shortest Mn-Mn separation stays nearly constant at 
around d2 ^2.85 A , the d\ distance shows remarkable 
agreement with Fig. [3j the ferromagnetic configuration 
shows stability over the helicoidal AFM one between 
50 K< T < 180 K where dx stays below 2.95-2.98 A. Al- 
though a direct extension of our K stability criteria to 
finite temperatures is ambitious, the temperature depen- 
dent magnetic behaviour of MnAso.88Po.i2 suggests that 
the criteria we have found can be applied more generally. 

The materials systems in this section exhibit a broad 
range oid2 values from 2.81 A in MnAso.88Po.i2 to 3.17 A 
in (Fei_xCoa^)MnP. Simultaneously, ^2 changes from be- 
ing the shortest to the second shortest Mn-Mn separa- 
tion. Despite this change, the variety of magnetic states 
found in these materials seem to be well described by crit- 
ical values of d\ only, as derived from our magnetic stabil- 
ity plot for MnP. In the following sections, we will further 
extend our analysis of this plot to XMnZ-type ternary 
alloys with the same orthorhombic structure. The im- 
portance of the d\ parameter in such materials will also 
become apparent. 



The first member of the series in which metamag- 
netism is observed is CoMnSi. A field- dependent meta- 
magnetic transition from a AFM helical state to a high 
magnetisation state was first reported in this composi- 
tion by Niziol and co- workers^^ . The transition temper- 
ature can be dependent on sample preparation, and was 
found to range from 207 K to 365 K. The peculiar and 
complex magnetic structure has been the subject of sev- 
eral invest igation^^^^. Early neutron diffraction studies 
found that a small magnetic moment is also developed on 
the Co atoms (~ 0.3/iB) in addition to the major mag- 
netic contribution from the Mn (~ 2.2/i^)^. The result- 
ing high saturation magnetisation (~120 emu/g at room 
temperature) and field-sensitive metamagnetic transition 
of the alloy makes it of interest as a potential room tem- 
perature magnetic refrigerant. We previously reported a 
magnetic field-induced isothermal entropy change (AiS) 
up tO'^ 7 JK'^kg"^ around room temperature for this 
allo}/^ and substituted variants^^ Recently, unusual 
negative thermal expansion along the a-axis of the lattice 
was also observed using high resolution neutron diffrac- 
tion^and in the absence of an applied magnetic field a 
crossover in di, and d2 is observed. In the above ex- 
amples, and in those to follow, the critical value of d\ 
parameter seems to predict well the change from FM 
to AFM order (around 2.95 A). This value, however, is 
probably slightly increased in the case of CoMnSi, and 
might be associated with the additional exchange com- 
petition among the Co-Co and Co-Mn pairs, compared 
to the binary MnP-type alloys. As stated above, Niziol 
et. al. identified up to 8 important magnetic interactions 
responsible for the overall magnetic properties^. 

The end member of the XMnSi series from FigjS] is 
NiMnSi, a ferromagnet with d^ =2.78 A and d2 =3.24 A 
at 80 K^^. The interpretation of the ferromagnetism of 
this ternary is apparent from the previous examples; a 
structure that accumulates the smaller Ni atom as a re- 
placement for Co brings the lattice into the zone where 
FM coupling of Mn-Mn atoms is more stable. In the fol- 
lowing section we show that metamagnetic behaviour in 
more general XMnZ-type Pnma structures can also be 
designed by the partial replacement of the p-block ele- 
ments. 



V. IMPLICATIONS FOR XMnZ-TYPE ALLOYS 
A. XMnSi 

There is a fairly large number of XMnSi-type al- 
loys that crystallize in the Pnma orthorhombic struc- 
tur d^^ * ^^^ ^^. However, many X elements are lanthanoids 
of large atomic size that result in considerable expan- 
sion of the lattice. The shortest and second shortest Mn- 
separations are plotted in Fig. [S] The high atomic ra- 
dius lanthanides cause a large expansion of the Pnma 
structure and increased Mn-Mn separation, stabilising 
the experimentally-found AFM states in accordance with 



B. NiMnGei-xSix 

Using the only ferromagnetic example from the pre- 
vious section, we now demonstrate that a stable AFM 
phase can be re-introduced even in the NiMnSi alloy by 
the increase of lattice dimensions with larger p-block el- 
ements. Just as with the pseudo-binary MnAsi_xPa; in 
subsection |IVB[ partial replacement of Si by Ge does 
not alter the overall valence electron number of the struc- 
ture. The NiMnGei_xSia; system was studied thoroughly 
by Bazela et. al.^^ and complex magnetic behaviour was 
observed. Three regions with different types of magnetic 
ordering were identified as a function Si concentration: 
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Figure 9: (Color online) Evolution of the two shortest Mn-Mn 
distances {di,d2) in pseudo-ternary NiMnGei-xS ix as a func- 
tion of composition at 80 K. Field-induced metamagnetism is 
reported within the range 0.25 < x < 0.55^. 



helicoidal AFM {0 < x < 0.25), AFM-FM type (NC) 
metamagnetic (0.3 ^ x < 0.55) and collinear FM with 
0.55 < X < 1. di and d2 are plotted as a function of 
Si concentration (top scale) and effective atomic radus, 
Reff (bottom scale) in Fig.|9] The lack of reported exper- 
imental error in the structural parameters derived from 
neutron diffraction data^^ only allow us to plot a linear 
fit line over the data points as a guide to the eye. The 
validity of the stability plot calculated for MnP in Fig. |3j 
however, is once again manifest for such pseudo-ternary 
NiMnGei_xSiaj compositions: a collinear FM state is sta- 
ble for di values below di ^2.95 A and gives way to AFM 
states for di ^2.95 A with increasing effective atomic size 
(Ge concentration). 



VI. CONCLUSIONS 

We have investigated the occurrence of AFM and FM 
states in Mn-based orthorhombic {Pnma^ 62) alloys. Us- 
ing DFT theoretical calculations based on the prototype 
MnP binary composition, a magnetic stability plot was 
constructed taking into account the two shortest Mn-Mn 
distances ((^1,(^2) • We have correlated literature obser- 
vations of different magnetic states in Mn-based Prima 
alloys with this plot and have found a remarkable agree- 
ment. As a result of isotropic expansion the FM ground 
state is no longer stable but instead AFM coupling of 
the spins on the Mn atoms is predicted when di ^2.95 A 
and d2 ^2.88 A, respectively. The di distance seems 
to play the dominant role in determining the magnetic 
state. In most of the samples the ^2 distance stays fairly 
constant as in the case of the (Co,Fe)MnP system or 



in the temperature dependence of the MnAso.88Po.12 

al- 
loy. In many other cases ^2 even shows the opposite 
tendency to di, as seen in NiMn(Si,Ge). The reason 
for these phenomena become clearer, if one examines the 
Pnma crystal structure in Fig. [T] The Mni-Mn2 atoms 
separated by the distance d2 form a continuous chain 
along the a- axis. This structural feature prefers spins on 
these Mn sites to be coupled antiferromagnetically by su- 
perexchange through occupied valence states along Mn- 
P-Mn. In fact, the magnetism along the chain probably 
favours AFM coupling even below the critical distance 
of d2 ~2.88 A but the theoretical calculations are based 
on an isotropically deformed lattice and the ratio of the 
a, b and c parameters for each e: value have been fixed. 
Further calculations with anisotropically strained lattice 
needs to be carried out in order to clarify this point. 
The Mn atoms separated by the nearest-neighbour dis- 
tance di (Mn3-Mn4), on the other hand, do not form 
a continuous chain in the lattice and so short-range di- 
rect exchange is the dominant interaction. Although the 
competition among direct and indirect exchanges that 
are also responsible for the observed non-collinearity in 
many of these alloys is foreseeable, the overall dominance 
of the direct exchange of Mn atoms along the di sepa- 
ration was not previously evident. The key role of this 
separation, resulting in FM coupling for di ^2.95 A and 
AFM coupling above di ^2.95 A and its control of the 
overall magnetism has been shown both theoretically and 
by comparison of many examples from the experimental 
literature. Interestingly, we find this critical di value, de- 
rived from binary MnP, to be influenced very little by the 
increased 4c site occupation present in Pnma ternaries. 

Based on these findings, it should be possible to design 
new alloys with lattice parameters tuned to a region of 
competition between AFM and FM states, enabling field- 
or temperature-driven metamagnetism. One example of 
such a material is FcxMnP with x ^ 0.6, where an Fe 
deficiency could reduce the lattice volume and bring di 
into the proximity of FM stability. Based on our findings 
we believe that there is still a large number of Mn-based 
ternary and pseudo-ternary alloys that can be created 
with metamagnetic transitions of potential interest for 
their magnetocaloric effect. 
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